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ABSTRACT 



Context. Small scale magnetic fields can be observed on the Sun in G-band filtergrams as MBPs (magnetic bright points) or identified 

in spectro-polarimetric measurements due to enhanced signals of Stokes profiles. These magnetic fields and their dynamics play a 

crucial role in understanding the coronal heating problem and also in surface dynamo models. MBPs can theoretically be described 

to evolve out of a patch of a solar photospheric magnetic field with values below the equipartition field strength by the so-called 

convective collapse model. After the collapse, the magnetic field of MBPs reaches a higher stable magnetic field level. 

Aims. The magnetic field strength distribution of small scale magnetic fields as seen by MBPs is infeiTed. Furthermore, we want to 

test the model of convective collapse and the theoretically predicted stable value of about 1300 G. 

Methods. We used four different data sets of high-resolution Hinode/SOT observations that were recorded simultaneously with the 

broadband filter device (G-band, Ca 11-H) and the spectro-polarimeter. To derive the magnetic field strength distribution of these small 

scale features, the spectropolarimeter (SP) data sets were treated by the Merlin inversion code. The four data sets comprise different 

solar surface types: active regions (a sunspot group and a region with pores), as well as quiet Sun. 

Results. In all four cases the obtained magnetic field strength distribution of MBPs is similar and shows peaks around 1300 G. This 

agrees well with the theoretical prediction of the convective collapse model. The resulting magnetic field strength distribution can be 

fitted in each case by a model consisting of log-noiTnal components. The important parameters, such as geometrical mean value and 

multiplicative standard deviation, are similar in all data sets, only the relative weighting of the components is different. 

Key words. Sun: magnetic topology. Sun: surface magnetism. Sun: atmosphere. Techniques: high angular resolution. Techniques: 
spectroscopic. Methods: observational 



1. Introduction 

The dynamics of the Sun's atmosphere are dominated by mag- 
netic fields. These fields span several or ders of magnitude in 
both field strength and in size (see e.g. IZwaanI [T987|). From 
extended features like sunspot groups (kG fields) to small ac- 
tive patches and pores, down to the smallest detectable el- 
ements, isolated single flux tubes. With newly installed ob- 
serving facilitie s (Hinode, NST, Su nrise; iKosugi etaD l2007l 
ICao et al.1 l2010l iBarthol et all l201lh and highly sophisticated 
methods like inversions of spectro-polarimetric data that are 
used more and more (e.g. the SIR code, the M erlin code; 
iRuiz Cobo & del Toro Iniestall992llLites et al.l2007h . flie detec- 
tion limit is constantly improving for weak and small magnetic 
fields in the solar atmosphere. Therefore, features smaller in size 
can be studied as can those weaker in magnetic field strength. 
This has led in recent years to a dramatic change in our under- 
standing of the quiet Sun's magnetic fields. The magnetic fields 
of the quiet Sun (which exclude active regions) comprise net- 
work and intranetwork fields. The network fields spatially cor- 
respond to the supergranular boundaries. Inside of these super- 
granu les, the fields are called intranetwork fields (see e.g. lZwaanI 
Il987h . 

While it was previously thought that the Sun's magnetic field 
in quiet regions is built up by strong vertical magnetic flux tubes 



in the intergranular lanes (between granules) and field-free re- 
gions in between (granules), recent studies have revealed that 
the field-free thought domains (intranetwork regions) are in fact 
filled up with ma gnetic fields with strengths up to hundreds 
of Gauss (see e.g. Orozco Suarez et al.ll2007l iLites et al]l2008l 



iDanilovic et a l."2010'). These intranetwork fields contain a large 
number of horizontal fields and can be observed above gran- 
ules. In many cases, horizontal intranetwork fields emerge co- 
spatially with granules a nd can be interpreted as low-lying loops 
(see e.g. Gomory et al. l2010h . However, there are also recent 
studies suggesting other interpretations of the new observations 
than strong hoiizontal magnetic fields (see e.g.: lAsensio RamosI 
120091 IStenfldll20To[|Borrero & Kobelll20TTl) . 

Such small-scale and weak magneti c fields can be ob- 
served in spectropolarimetric data (e.g . Bello Gonzalez et al] 



. e.g. mello 

120081 iNagata et al.ll2008L IViticchie et al.ll2009h or in magnetic- 
field sensitive filtergrams, such as the G-band (e.g.i Bergeret al.l 
1998, 2004). In such filtergrams, the magnetic field concentra- 
tions ar e identified as so-called magnetic bright points (MBPs; 
see e.g. |K eneiill992l lYi & Engvoldll 19931 iBerger & Titlell200lL 
ISanchez A lmeida et al. 2004 lUtz et al.1 l2009al) . The G-band 
(centred on 430.5 nm) is primarily used for such investigations 
because it gives a better contrast between MBPs and the sur- 
rounding granulation than a continuum image would. The in- 
creased contrast is due to a partial evacuation of the flux tube. 
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Fig. 1. Magnetic field strength map of the fully inverted SP data of data set I, together with the identified MBPs (indicated by 
crosses). The white circle outlines an exclusion zone around the sunspot; i.e. identified MBP features within this zone are not taken 
into account in the analysis. 



This leads to lower opacity, hence a deeper formation height in 
the atmosphere. As we can see deeper and therefore in a hotter 
atmospheric layer, the opacity in the G-band is further decreased 
owing to an increased dissociation rate of the contributing CH 
molecules. This gives rise to a higher contrast of MBPs in the 
G-band compared to the continuum. Fo r mor e details about 
these processes, see ISteiner et aLl (1200 lb and ISchiissler et al.l 
(l2003h . The relationship between the contrast and the magnetic 
field strength was confirmed by simulations like the one by 
IShelvagetal.l(l2004l). 

In the 1970s, ISpruitJ (Il979h developed a theoretical model 
for the formation process of MBPs, the so-called convective 
collapse model. The model suggests that when magnetic fields 
exceed a critical magnetic field strength (equipartition field 
strengthQ), the plasma within the magnetic field cools down to 
form down drafts within the field, which finally leads to a col- 
lapse of the magnetic field, resulting in a small scale vertical flux 
tube. This flux tube is then visible as MBP. The model suggests 
an equipartition magnetic field strength after the collapse of the 
magnetic field to a flux tube of roughly 1300 G. 

In this paper, our aim is to study the magnetic field strength 
distribution of small scale magnetic fie lds observed as MBPs. A 
similar investigation was conducted bv lBovelet & W iehr (20081). 
These authors identified and analysed MBPs in G-band filter- 
gram data to obtain the magnetic flux contained in these small 
scale kG field strength exhibiting features of the solar photo- 
sphere. For estimating of the magneti c field strength, they used 
a calibration curve of IShelvag et al.l (|2004|) . Nevertheless the 
direct estimation of the magnetic field strength distribution of 
MBPs on a larger scale and for different solar surface regions 
was not done beforq^ due to the lack of high-resolution co- 



' The field strength at which the magnetic pressure has reached a 
level sufficient to balance the kinetic pressure exerted on the mangetic 
field; so the magnetic field cannot get further compressed by the sur- 
rounding plasma flows. 

^ There is a very interesting study bv lBeck et al.l(l2007h of MBP prop- 
erties in the sunspot moat. Among the studied parameters is the mag- 
netic field strength and flux. 
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Fig. 2. The brightness excess distribution of MBPs identified in 
the four data sets. The brightness excess is defined as the bright- 
ness at the brightness barycentre position of an MBP divided by 
the mean brightness of the FOV (both taken from the G-band 
images). 



temporal data sets comprising spectro-polarimetric data (suit- 
able for inversions to gain magnetic field strength maps) and 
filtergram data (especially G-band data for the identification 
of MBPs). The Hinode satellite deployed in 2006 provides the 
right data for such an investigation since it combines a filter- 
gram imager (G-band observations) with a high-resolution spec- 
tropolarimetric device (SP - spectro-polarimeter). Details about 
the Hinode mission and its SOT (Solar Opt ical Telescope) in- 
strume nt can be found in lKosugi et al.l (l2007h and lTsuneta et al.l 
(I2008h . 



2. Data sets 

We studied four different data sets comprising fully inverted SP 
spectropolarimeter data and G-band filtergrams. The SP inver- 
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Table 1. Fit parameters (cr and fi) and weighting coefficients (wi and W2) derived for tfie magnetic field strength distribution of data 
set I (sunspot group; see Figs.[T]|3]l. The reduced ;t'^ value (describing the goodness of the fit) has a value of 3.9. 



description 

log-normal distribution component I 
log-normal distribution component II 



weighting coefficient relative weighting coefficient mean value (/j) standard deviation (cr) 

[%] [G] [G] 



72000 ± 2000 


66 + 2 


5.62 + 0.03 


0.80 ± 0.02 


37000 ± 1500 


34 ±2 


7.13 ±0.01 


0.17 ±0.01 
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FOV the variation in the inclination angles is within 4. 1° to 17.4° 
(cos 6 value of 0.99 to 0.95), which is low enough that no spe- 
cial care has to be taken for projection effects. The inverted SP 
data (magnetogram) has a spatial sampling of 0.30 arcsec/pixel 
in the x-direction and of 0.32 arcsec/pixel in the y-direction. The 
coiTesponding filtergram data set had the following characteris- 
tics. The temporal resolution for the filtergram images is about 
seven minutes. The spatial sampling was reduced by onboard 
binning to 0.108 arcsec for the G-band and CaiiH filtergrams. 
In addition, Fe line filtergrams of the narrow band imager (NFI) 
instrument are available. These data have a spatial sampling of 
0.16 arcsec/pixel. The filtergram time series covers the peiiod 
from 15:32 UT to 18:58 UT with a different FOV of 221 arcsec 
by 110.5 arcsec. 



Fig. 3. The magnetic field strength distribution of MBP features 
gained by data set I (sunspot). The peak around 1 300 G related to 
collapsed magnetic fields was fitted by a log-normal component 
yielding a geometric mean yu = (1250 ■ |h- 1.18)G. The reduced 
X^ for the fit is 3.9. In total, 2392 MBP features were identified 
and used in the analysis which covered about 2 % of the FOV. 

sions were carried out with the Merlin code|_| (courtesey to Bruce 
Lites and co-workers). The data sets contain different typical so- 
lar features, suchas a sunspot group, a magnetically active region 
(pore), and quiet Sun. These data sets were chosen to a) cover 
different features of the Sun and b) feature nearly co-temporal 
filtergram and spectropolarimetric data (except data set I) with a 
high temporal cadence. 

Data Set I (sunspot group): The data stem from 2006 
December 1 1 and show a sunspot group (see Fig.[T]Q A few days 
later this active region gained a reputation bec ause it hosted the 
first X-cl ass flare observed by Hinode (see e.g. iKubo et al]|2007l 
iJing et al. 2008. Asai et al. 2008) . In addition to flie SP data, data 
of the BFI (broad band filter imager) and NFI (naiTow band fil- 
ter imager) device are available (G-band, CaiiH, Fe line). The 
main features of the SP data are as follows. The data were ob- 
tained by scanning the surface in the fast-map mode and cover 
the period from 17:00 UT to 18:03 UT. The totally scanned field 
of view (FOV) has a size of 295.2 arcsec by 162.3 arcsec. The 
centre of the FOV pointed to solar coordinates of x = 152.4 
arcsec, y = -96.1 arcsec. This position on the Sun corresponds 
to a heliocentric angle (0, inclination of vertical flux tubes to the 
line of sight) of 10.6° equaling a cos lvalue of 0.98. For the total 



■^ Detailed information about the Merlin code can be 
found in the internet under the following web URL: 
[http ://www.csac.hao.ucar.edu/csac/nextGeneration.jsp 

^Some pixels in the centre of the sunspot were not correctly inverted. 
This is most probably due to the strong magnetic fields in the centre of 
the sunspot. For our study this is of no concern as we are interested 
in the magnetic field strength of the MBP features which are located 
outside of the sunspot. 



Data Set II (active region; pores); Data set II covers a small 
active region outside of the disc centre (shifted to the noith-east; 
X = -330 arcsec and y - 390 arcsec). In this case the helio- 
centric angle has on average a value of 32.7° with a variation 
within 30.2° to 35.3° corresponding to a mean cos 0of 0.84±0.2. 
The exposures were taken over a time interval of about 18 min- 
utes with a one-minute cadence for the filtergram exposures and 
a two-minutes cadence for the SP scans. The time series staits 
at 9:00 UT on 2009 June 2 and lasts until 9:20 UT. The data 
comprise Hinode SOT/SP data fully inverted by the Merlin code 
(courtesy to B. Lites et al.), G-band, Ca 11 H observations and Fe 
I line exposures of the NFI filter The FOV of the filtergram ex- 
posures is about 80 by 90 arcsec^. The SP scans show an FOV 
of 9.45 arcsec by 81.15 arcsec and compiise ten fast map scans 
of the region. The spatial sampling is 0.295 arcsec/pixel. The 
Fe I line filtergrams are available in all four Stokes parameters 
(/, Q, U, V) with a reduced FOV (compared to the BFI filter- 
grams) of 12.8 arcsec by 81.92 arcsec. The spatial sampling rate 
is higher than for data set I and amounts to 0.16 arcsec/pixel. 
The time series for the filtergram exposures covers a longer pe- 
riod, from 8:25 UT to 9:55 UT. In total, the filtergram data set 
comprise 90 exposures in each wavelength. 



Data Set III (quiet Sun); Data set III was recorded on 2007 June 
2. It covers a quiet Sun region near disc centre (about 200 arcsec 
south of the disc centre). This gives a heliocentric angle for the 
FOV, hence also for the flux tubes in a range of 8.8° to 15.6° with 
a mean value of about 12.2°. These values correspond to a mean 
cos value of 0.98. For the total FOV the values range between 
0.96 and 0.99. The data comprises SOT G-band, CaiiH, and SP 
data. The SP data were inverted by the Merlin code. The tempo- 
ral coverage is about 2 minutes for the magnetograms and about 
30 seconds for the filtergrams. In total this gives 26 fast-map SP 
scans of the region with a spatial sampling of 0.295 arcsec/pixel 
in the time span from 11:48 UT to 12:40 UT The FOV of the 
SP data was 8.86 arcsec by 162.3 arcsec. Filtergram exposures 
were taken with a spatial sampling of 0.108 arcsec/pixel starting 
at the same time but lasting until 12:41. This gives 106 images 
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Fig. 4. Top: typical G-band filtergram of data set II (active re- 
gion). The region scanned by the SP instrument is marked by a 
rectangle. Bottom: available data and identified MBPs (marked 
by white crosses). From left to right: whitelight image (obtained 
by inversion; SP data), Ca II-H filtergram (chromosphere), mag- 
netic field strength map (gained by inversion; SP data), and 
finally a G-band filtergram. The range of the magnetic field 
strengths for this data is given as 14 to 2580 G. 



for each of the two filters used. The FOV of the filtergrams is 
55.8 arcsec by 1 1 1.6 arcsec. 

Data Set IV (quiet Sun, HR data): The major diff'erence of 
this data set to the previous ones is the higher spatial sam- 
pling. The instrument always obtains the same/or nearly the 
same diffraction limited data sets. Nevertheless, the resolution 
quality can be slightly c hanged by the spatial sampling rate 
(see e.g. lUtz et al]|2009bh . While in the case of the other data 
sets, the data were binned two by two, data set IV comprises 
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Fig. 5. The magnetic field strength distribution obtained for data 
set II (active Sun; pores). The fit was created by using two 
log-normal components. The first component assigned to MBPs 
yielded a geometric mean value of (1260 ■ | H- 1.25) G, while the 
second component is used for the background, which is rather 
weak. The reduced ;t'^ for the fit is 3.0. For this analysis, 2624 
MBP features were taken into account and occupied a fraction 
of the FOV of about 6%. 



unbinned data. Therefore, the spatial sampling is doubled with 
0.054 arcsec/pixel for filtergram images and 0. 16 arcsec/pixel for 
the spectropolarimeter This is the best spatial sampling Hinode 
can provide for the SOT instrument. The data were obtained on 
2007 January 20 and were recorded between 6:54 UT and 7:58 
UT, so they span roughly one hour of solar quiet Sun evolution. 
The G-band filtergram data were taken with a temporal reso- 
lution of 30 seconds, whereas the SP data have a cadence of 
roughly two minutes, so about four filtergram exposures belong 
to one SP slit scan. The FOV of the filtergram images is 400 
pixels by 1024 pixels or 21.6 arcsec by 55.3 arcsec, respectively. 
The FOV of the SP data is 25 pixels by 5 12 pixels, corresponding 
to 3.7 arcsec by 81.9 arcsec. The SP instrument has been oper- 
ated in normal mode compared to fast-map mode in the other 
data sets, which yields a higher spatial sampling of the data. The 
selected FOV (same FOV for BFI and SP) is about 3.7 arcsec by 
55.1 arcsec. The SP device took in total 29 scans, whereas the 
filtergram imager took 128 exposures. The centre of the FOV 
had solar coordinates of .x = -4.5 arcsec shifting during the time 
series to 5.6 arcsec and a y =2.9 arcsec. We can therefore as- 
sume, in very good approximation, that these data were all taken 
at the solar disc centre (heliocentric angle of 0° and cos with a 
value of 1). The inclination of the flux tubes to the angle of sight 
therefore has (if at all) a negligible influence on the results. 



3. Analysis 

To analyse the magnetic field strength distribution of MBPs it 
is necessary to identify the MBPs correctly and then to deter- 
mine their corresponding magnetic field strength. After obtain- 
ing beneficial datcQ (see Sect. 2) from the Hinode data base, the 



' These data sets must, on the one hand, comprise at least longitu- 
dinal magnetic field strength maps (coming e.g. from inversions; in the 
following, we refer to the magnetic field strength map somewhat loosely 
as a magnetogram) and on the other hand filtergram data of at least the 
G-band filter. Furthermore the temporal and spatial resolution should 
be high, and the magnetogram data to the filtergram data has to be as 
co-temporal as possible. 
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Table 2. Fit parameters for the MBP field strength distribution of data set II (active region; Fig.|5]l. The fit gives a reduced ;\f^ value 
of 3.0. 



description 

upper panel; log-normal component I 
upper panel; log-normal component II 



weighting coefficient relative weighting coefficient mean value (jj) standard deviation {cr) 

[%] [G] [G] 



40000 ± 6000 


31±5 


6.6 ±0.1 


0.59 ± 0.05 


87000 ± 6000 


69 ±5 


7.14 ±0.01 


0.22 ± 0.01 



Table 3. Fit parameters {cr and fi) and weighting coefficients {w[ and W2) derived for the magnetic field strength distribution of data 
set III (quiet Sun; see Figs. |6l|7]l. The reduced ;i'^ goodness of fit parameter yielded a value of 8.3. 



description 

log normal distribution component I 
log normal distribution component II 



weighting coefficient relative weighting coefficient mean value (ju) standard deviation (cr) 

[%] [G] [G] 



116000 ±2000 


81 ± 1 


5.69 ±0.01 


0.82 ±0.01 


28000 ± 1000 


19 ±1 


7.208 ± 0.003 


0.085 ± 0.002 
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Fig. 7. Magnetic field strength distribution of the identified MBP 
features in a quiet solar FOV. The hump related to MBPs is seen 
at (1350 -K 1.09) G. The quahty of the fit (solid fine) is 8.3 
(reduced x^)- The distribution was created by analysing 6253 
identified MBP features coveting about 1 % of the FOV. 



Fig. 6. From left to right: Continuum image stemming from SP 
data, Call-H (chromosphere), magnetic field strength map (in- 
version of SP data done by MERLIN code), G-band (photo- 
sphere). The crosses indicate identified MBPs. It can be clearly 
seen that no significant magnetic field activity is within the FOV 
(pores, sunspots). The magnetic field strengths in the data set 
vary between 3 and 1510 G. 



next step is a careful alignment of the filtergram and magne- 
togram data sets. The G-band filtergrams have a different spatial 
and temporal sampling than the fully inverted level-2 spectro- 
polarimetric data. We overcame the different spatial sampling by 
stretching the magnetogram data sets to the corresponding spa- 
tial sampling of the filtergrams. This can be done in a first step by 
roughly expanding the magnetogram data set by the correspond- 
ing spatial sampling ratio. Since MBPs are very small-scale fea- 
tures, a rough alignment of + 1 arcsec would be insufficient. For a 
better alignment one can use the continuum intensity within the 
IeveI-2 SP data sets. These continuum images show a granulation 
pattern similar to the G-band. By calculating cross-correlation 
coefficients, the expanded image can be aligned to the magne- 
togram on pixel and even subpixel scales. Finally, the data are 



cut to the same FOV. For details about high precision alignment 
and a coiTesponding algorithm for Hinode/SOT data, we refer to 
Kuehner et al. (2010). The temporal alignment was done by tak- 
ing the temporally closest filtergram to an SP scan. This gives 
temporal alignments with differences smaller than 0.5 minutes 
for data sets II to IV. In the case of data set I, we made just a 
spatial alignment of one filtergram to the complete SP scan. For 
problems evolving out of this we refer to the discussion section. 
After careful alignment the next step is to identify the in- 
teresting MBP features within the G-band data. This was done 
by the algorithm of Utz et al. (2009a). The algorithm takes in a 
first step the G-band images and segments them into individual 
features. The idea of the segmentation is to follow the biight- 
ness contour of the features within an image from the brightest 
pixels down to the faintest. After applying of the segmentation 
the interesting MBP features are identified. The version applied 
here differs to the version of the cited paper (where the identi- 
fication was based on the size of the features) in the sense that 
for the identification the local biightness gradient is considered. 
This means that for each segment the brightness ratio between 
the brightest pixel of a segment and an average brightness in 
a prescribed vicinity is calculated. If this ratio is high enough, 
the feature is identified as an MBP. If the ratio is close to 1, the 
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Table 4. Fit parameters (cr and //) and weighting coefficients (wi and W2) derived for the magnetic field strength distribution of data 
set IV (quiet Sun, high-resolution data set; see Fig.|8]l. The measured goodness of fit in this case was 4.4. 



description 



weigiiting coefficient relative weighting coefficient mean value (/j) standard deviation (cr) 

[%] [G] [G] 



log-normal distribution component I 78000 ± 1000 

log-normal distribution component II 6300 ± 500 



93 + 2 
7± 1 



5.98: 
7.05: 



0.01 
0.01 



0.52: 
0.10: 



0.01 
0.01 



segment belongs to a granule. The global brightness exces^ dis- 
tribution of MBPs is shown in Fig. |2l The rest of the analysis 
is quite straightforward. After identifying the correct features, 
the barycentre position of the MBPs is calculated. The coordi- 
nates of the barycentre are then used to extract the coiTespond- 
ing magnetic field strength in the aligned magnetograms. The 
obtained magnetic field strengths are then represented in a his- 
tograrrQ, and the resulting distribution was finally fitted by a two- 
component log-normal distribution. The fit equation therefore is: 



Magnetic field strenght distribution 



f{o-,fi,x) = Wigio-util,x) + W2gia-2,IU2,x), 



(1) 



where gicr,fi, x) stands for the so-called log-normal distribution, 
W] and W2 specify weighting coefficients. These weighting coef- 
ficients give the relative contribution of each component to the 
total function. The log-normal distribution g{cr, fj., x) is defined 

as 

' log X — fl\ 



g{o-,n,x) 



1 1 

xcrSln \ 2 



2\ 



(2) 



The parameters are cr, a multiplicative standard deviation, 
and ji the geometric mean value. These parameters are some- 
what different from the generally known yU and cr parameters 
used in normal distributions as we empl oy log-normal distri bu- 
tions (more information can be found in iLimpert et al.ll2008h . 

The extraction of magnetic field information from magne- 
tograms could lead to some errors for strongly inclined flux 
tubes. In our case the approach is appropriate because a) the data 
sets are close to the solar disc centre, hence producing low incli- 
nation angles, and b) MBPs are in general connected with strong 
vertical flux tubes. 



4. Results 

Data Set I (sunspot group): The obtained magnetic field 
strength distribution for MBP features is shown in Fig. |3] In the 
figure caption of this figure and the following figures, the symbol 
■|-^ means that the one cr boundary can be gained by multiplying 
(higher boundary) and dividing (lower boundary) the geometric 
mean value by the multiplicative deviation. We use log-normal 
distributions and not the commonly used normal distributions 
for which the boundaries can be calculated by adding and sub- 
tracting the standard deviation. The fit values we derived are 



* Global means in this case the excess between an MBP feature to 
the mean G-band brightness of the total FOV. For the identification only 
local brightness excesses were considered. 

' The binsizes of the histograms were calculated by dividing the 
range of magnetic field strengths by an adequate number of represent- 
ing classes. The number of classes can be estimated by calculating the 
square root of the number of data points, i.e. that the number of bins 
(and also the width of the bins) can vary from data set to data set de- 
pending on the total number of identified feature s. For the problem of 
correc tly binning histograms see the treatise in ISalgado-Uearte et al.l 
( I2OOOI) comprising examples and clarifying the impact of changing bin- 
widths. 




1000 1500 

Field strength [G] 



Fig. 8. Magnetic field strength distiibution of the identified MBP 
features in the quiet solar FOV of the high spatial resolution data 
set (data set IV). The reduced;^'^ for the fit is 4.4. The component 
related to MBPs is fitted at (1110 -In- 1.14 ) G. For this data 
set, 2614 identified MBP features were analysed, which covered 
about 1 .4 % of the total FOV. 



listed in Table [T] The two fitted components correspond to one 
log-normal component that fits the magnetic field strength back- 
ground signal and a second one that fits flux tubes with higher 
magnetic field strengths. This might be MBPs in their collapsed 
state; post-collapse MBPs. The background distribution can be 
explained by considering of MBPs before they undergo the col- 
lapse, during the collapse, and probably during the dissolution. 
Also a fraction of wrongly identified, non-MBP features might 
be taken into account. Finally the shape of this component is re- 
lated to the dynamo process itself, creating magnetic fields up 
to the equipartition field strength, and to the way instrumental 
noise contributes via the inversion tools to the magnetograms. 
The more interesting component belongs to the collapsed flux 
tubes and is centred at 1250 G with a multiplicative standard 
deviation of 1.18 and a relative contribution of 34%. This com- 
ponent can be attributed to the field strength of MBPs after the 
convective collaps e. 1350 G is p redicted by the convective col- 
lapse model e.g. in lSpruitl(II979l) . 

In the case of this data set one should have in mind that the 
SP instrument scans over the FOV. Such a scan needs about one 
hour for the actual size of the FOV. Therefore the left-hand side 
of the image shows the Sun's magnetic field at about 17:00 UT, 
whereas the right-hand side of the image shows the Sun's mag- 
netic field at about 18:00 UT. We keep in mind that the evolution 
of small-scale fields happens in th e range of minutes (see e.g.: 
Ide Wiin et aI.ll2005llUtz et al.l2010h . Therefore it is clear that the 
correlation between a single G-band filtergram and the magne- 
togram can only be statistical, i.e. that only a fraction of MBPs 
will be measured in the collapsed state with the correct mag- 
netic field (giving rise to a relative contribution of 34%), while 
for the others the magnetic field before or after their appearance 
will be measured (giving rise to a broad background log-normal 
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Table 5. Parameters for the log-normal component of the MBP magnetic field strength fit for all four data sets of our study. The 
mean value in a log-normal distribution is the geometrical mean and that the standard deviation is not addivitive but multiplicative. 



description 


mean value fit (/i) 


mean value phys. 
[G] 


standard deviation fit (cr) 


standard deviation as 
multiplicative factor 


data set I (large FOV with sunspot) 


7.13 ±0.01 


1250 ±10 


0.17 ±0.01 


1.18 ±0.02 


data set II (active Sun with pores) 


7.14 ±0.01 


1260 ±10 


0.22 ±0.01 


1.25 ±0.01 


data set III (quiet Sun) 


7.208 ± 0.003 


1350 ±10 


0.085 ± 0.002 


1.08 ±0.01 


data set IV (quiet Sun; HR data) 


7.05 ± O.OI 


I150±I0 


0.10 ±0.01 


1.10 ±0.01 



distribution component of 66%). For more details about the im- 
plications we refer to the discussion section. 



nent shows values of ju = 5.98 and cr = 0.52. For more details 
about the obtained parameters we refer to Table 



Data Set II (active region; pores): The main difficulty arises 
in the temporal and spatial alignment. The temporally closest 
subfields of the filtergrams corresponding to the scanned FOV 
of the SP instrument were selected and aligned (accuracy better 
than one pixel and two minutes, respectively). For an illustrative 
example see Fig.|4] Figure |5] shows the obtained magnetic field 
strength distribution for this data set. The fit follows again the 
two-component log normal model. The geometrical mean for the 
MBP related component has a value of (1260 ± 10) G and a 
multiplicative standard deviation of 1.25. All obtained fit values 
are given in Table |2] Both components have similar slopes and 
similar mean values compared to the case of data set I (sunspot) 
and yet the relative contributions of background and collapsed 
flux tubes have significantly changed from 66% to 3 1 % and 34% 
to 69%. 



Data Set III (quiet Sun): Figure |6] displays the different filter- 
gram and magnetogram data of the used data set, and Fig. Qgives 
the measured magnetic field strength distribution for this case. 
The fit shown in Fig. Q was created by the same fit function as 
for data set I, i.e. two log-normal components. The main finding 
compared to the other data sets is that the relative importance of 
the convective collapse magnetic field strength peak about 1 350 
G is relatively weaker than in the other data sets. For data set 
III the relative weighting value derived is only 19% compared to 
34% for data set I and 69% for data set II. The peak is centred 
at 1350 G, which is a slightly higher value compared to the one 
obtained for data set I. The variance of the log-normal compo- 
nent is 1.09 (multiplicative vaiiance of 9%), which is by a factor 
of 2 smaller compared to data set I for which we derived a value 
of 1.18 (multiplicative variance of 18%). For more details about 
the obtained fit values we refer to Table |3] 



Data Set IV (quiet Sun; liigli-resolution data): The analysis of 
this data set differs from the analysis of the other data sets by 
the exclusion of MBP features with sizes smaller than at least 
four pixels. This is justified for two reasons: a) the smallest sizes 
are more likely to contain noise and/or are produced by noise 
(esepecially due to the higher spatial sampling) and b) a four- 
pixel size represents the same minimum size detectable in the 
other data sets (binning 2 by 2), hence an exclusion of smaller 
sizes improves the comparability between the different data sets. 
Figure[8]shows the obtained magnetic field strength distribution 
for this data set and the results of the applied fitting model. As in 
the other data sets a significant log-normal component for col- 
lapsed flux tubes (MBPs) shows up at a value of (11 50 + 10) G 
for the high-resolution data. The multiplicative variance factor 
for this component is 1.10. The background log-normal compo- 



Summary: In all four data sets we were able to fit the ob- 
tained magnetic field strength distribution of MBPs with a model 
consisting of two log-normal distiibutions. One component was 
needed for the background (non-collapsed flux tubes) and the 
second component for the collapsed state of flux tubes. This type 
of distribution is known to be formed in agglomeration processes 
like the growing of molecules on substrates, e.g. size distribu- 
tions (see Limper t et al.l l2008h . In general, a random variable 
or measured quality shows this kind of distiibution if the un- 
derlying process, which leads to the formation of the value, is 
created by multiplicative processes rather than by additive pro- 
cesses. Therefore the obtained fit parameters and shapes of the 
distributions could yield a statistically significant insight into the 
magnetic field generation (dynamo) process. 

For detailed values of the fits see Table |5] The values reveal 
that the spread of the geometric mean for the MBPs ranges from 
1110 G to 1350 G (data set IV and data set III, respectively). 
This can be solved by the following idea, on which we elabo- 
rate more in the next section, that we will indeed not need one 
component to fit the collapsed state but rather two components. 
A weak component may be related to intranetwork fields peak- 
ing around 1100 G and a stronger component probably related 
to network fields peaking around 1300 G. With this in mind, the 
explanation for the range of results may be the different contri- 
butions of the network/intranetwork components in the four data 
sets. 

Neveitheless, at this point in our analysis we can conclude 
that the MBP field strength distribution follows log-normal dis- 
tributions with maxima between 1 100 G and 1400 G. These val- 
ues are wit hin theoretica l predictions of the convective collapse 
model (e.g. lSpruillll976h . 

5. Discussion 

We note that we are discussing magnetic field strengths in this 
paper, while other authors state that observations can only lead 
to the magnetic flux contained within a pixel. This means that the 
true magnetic field strength would depend on the filling factor 
We implicitly assume a filling factor for MBPs of one, i.e. that 
the pixel is filled completely with the magnetic field. This may 
be justified by the fact that MBPs stretch over several pixels and 
by the congruence between theoretical predictions and these ob- 
servations (the magnetic field strength value after the convective 
collapse of about 1300 G). On the other hand, it may explain 
how it is possible that many of the MBPs are measured with 
lower magnetic field strengths except of the active region data 
set, implying that the filling factor is thus in reality smaller than 
one and that the true structures are smaller and/or fragmented 
and yet not resolved. Nevertheless, these are speculations, so we 
will have to wait for future instruments with higher resolutions 
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Fig. 9. Top panel: simplistic fit model (two log-normal compo- 
nents) applied to data set I (red, short dashed line) compared to a 
medium complex model consisting of 3 log-normal components 
(purple, long dashed line) and a sophisticated (4 log-normal 
components) model blue, short/long dashed line). Second panel: 
illustration of the components of the simplistic model. Third 
panel: components of the medium complex model. Fourth panel: 
components of the sophisticated fitting model, which assumes 
two strong magnetic field components around 1000 G and 1300 
G, respectively. The reduced ;^^ value for the simplistic model is 
3.9 compared to 1.3 for the medium complex model and 0.8 for 
the most sophisticated one. 



to bring the true physical nature to light. A possible extension 
of this study can hopefully be do ne in the future with the IM aX 
instrument (for more details see Martinez Pillet et al. '201 1") on- 
board the balloon-borne Sunrise mission. Lagg et al. ( 2010 ) has 
shown in a recent study that quiet Sun magnetic fields can be 
resolved with this instrument. 



Fig. 10. Same as in Fig.|9]but for data set III. It can be seen that 
the components are in position (but not strength; weighting) sim- 
ilar to the ones in Fig.|9] The goodness-of-fit values (reduced;t'^) 
are 8.3 for the simplistic model (2 components; second panel) 
compared to 3.2 for the medium complex model (3 components; 
third panel) and 1.8 for the more sophisticated one (4 compo- 
nents; final panel). 

In all cases we were able to fit the background magnetic 
field distribution with a log-normal distribution. This distribu- 
tion consists of still uncollapsed flux tubes, intermediate states 
in the evolution of the flux tubes, and a low fraction of wrongly 
identified features (identification errors). Log-normal distribu- 
tions occur naturally by physical mechanisms involved in gener- 
ating of the magnetic field. On the other hand, it is thinkable that 
the background distribution is also generated to a certain extent 
by measurement noiso Especially the position and height of 



** This noise is composed of the photon noise (the higher the spectral 
resolution the lower the photon yield) and the noise induced by the ex- 
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Fig. 11. Same as in Fig.|9]but for data set IV. The goodness-of- 
fit value x^ is 4.4 for the simplistic model and improves over 
2.1 for the medium complex one to a final value of 1.9 for the 
sophisticated model. 



Fig. 12. Same as in Fig. |9]but for data set II. The simplistic fit 
gives a rather good agreement with the measurements (reduced 
X^ of 3.0); the medium complex model and the sophisticated 
model fit the data better with values of 1.5 and 1.1, respectively. 



the peak of the background distribution, which lies in our case 
between 100 G and 200 G, is quite often repo rted as a nois e 
arte fact, and it shifts with image quality (see e.g.lS tenfloll2010h . 
ISten flo (2010) states, for example, that if one corrected the influ- 
ence of noise on the measurements of quiet Sun magnetic fields, 
there would be no peak at all, and the probability distribution 
should rise very strongly for field strengths approaching zero. 
Furthermore, it cannot be excluded that this peak at about 190 
G is caused by the interaction of the detectability (image noise) 
of magnetic fields and the inversion codes, which transform the 
polarisation signals in magnetic field vectors. 



posure time: the longer the exposure time the better the photon statistics 
but the worse the contrast of the image due to the temporal averaging 
over the dynamics of the features contained within the FOV. 



Data alignment: The alignment of filtergram data with spectro- 
polarimetric data inhibits an intrinsic temporal alignment prob- 
lem. This is due to the design of the instrument as slit spectro- 
graph. Therefore only data from a tiny portion of the Sun (one 
slit width) can be gained instantly. To obtain a suitable FOV, the 
slit must be moved over the surface to be mapped. This scanning 
process takes some time depending on the size of the FOV. This 
problem can be circumvented by other instrumental desig ns (see 
e.g. the Sunrise mission with the IMaX instrument; Barth ol et al.l 
1201 ll iMartinez PiUet eFani201 ll) . which naturally inhibits other 
drawbacks (such as lower spectral resolution). Magnetograms 
taken by a slit spectrograph (like Hinode/SOT/SP) therefore 
show a temporal evolution of the Sun's surface. 
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Table 6. Parameters for the medium complex fitting model for the magnetic field strength distribution fits as shown in Figs.l9lfT0l 
fm and[T2] "Data set" specifies the used data set; w, cr, and/v specify the log-normal components; and;\f^ represents the goodness- 
of-fit value. The indexing number indicates the corresponding component. The parameters of this table should be compared with 
the parameters of the simplistic model (as stated in Tables[T]|2l[3]|ll or the sophisticated model as summarised in Table |7] 



data set 


(L>1 


^1 


O"! 


C02 


A'2 


0"2 


OJi 


/^3 


0"3 


T 

X' 


sunspot 

active 

quiet 

quiet HR 


54 ±2 
10 ±2 
51±2 
64 ±3 


5.46 ± 0.04 

6.0 ± 0.2 

5.25 ± 0.04 

5.80 ± 0.03 


0.80 ± 0.03 

0.8 ±0.1 
0.67 ± 0.02 
0.51 ±0.02 


25 ±2 
60 ±3 
18 ± 1 

7± 1 


7.20 ± 0.03 

7.18 ±0.01 

7.212 ± 0.003 

7.10 ±0.01 


0.13 ±0.01 

0.21 ±0.01 

0.083 ± 0.003 

0.09 ± 0.01 


21 ±2 
30 ±3 
31±2 
29 ±3 


6.72 ± 0.03 
6.74 ± 0.04 
6.49 ± 0.02 
6.47 ± 0.03 


0.31 ±0.03 
0.37 ± 0.03 
0.39 ± 0.02 
0.30 ± 0.02 


1.3 
1.5 

3.2 
2.1 


averaged 


- 


5.6 ±0.3 


0.7 ±0.1 


- 


7.17 ±0.05 


0.12 ±0.06 


- 


6.61 ±0.15 


0.34 ± 0.05 


- 



Table 7. Parameters for the sophisticated magnetic field strength distribution fits as shown in Fig. |9] [JO] [TT] and [12] "Data set" 
specifies the used data set, x^ the goodness of fit, co, cr and jj specifies the log-normal components. This table should be compared 
with the previous tables of the simplistic model ([Tl|2]|3]|4]l or the intermediate model (see Table|6}. 



data set 


x' 


component 


CO 


^ 


cr 


component 


tLl 


^^ 


cr 


sunspot 


0.8 


1 

3 


48 ±5 
18±7 


5.39 ± 0.09 
6.5 ±0.1 


0.77 ± 0.04 
0.4 ±0.1 


2 
4 


20 ±5 
14 ±7 


7.24 ± 0.01 
6.97 ± 0.07 


0.11 ±0.01 
0.18 ±0.05 


active 


1.1 


1 
3 


11 ±20 
21 ±9 


6± 1 
6.6 ± 0.3 


0.9 ± 0.3 
0.4 ± 0.2 


2 
4 


10 ±6 

58 ± 14 


7.31 ±0.02 
7.12 ±0.03 


0.10 ±0.03 
0.23 ± 0.02 


quiet 


1.8 


1 

3 


43 ±2 
36 ±1 


5.12 ±0.04 
6.33 ± 0.03 


0.62 ± 0.02 
0.42 ± 0.03 


2 
4 


16 ±2 
8±3 


7.23 ±0.01 
7.00 ± 0.04 


0.073 ± 0.004 
0.14 ±0.04 


quiet HR 


1.9 


1 

3 


64 ±3 

27 ±3 


5.80 ± 0.03 
6.44 ± 0.03 


0.51 ±0.02 
0.29 ± 0.03 


2 
4 


4± 1 

5±2 


7.16 ±0.01 
6.99 ± 0.03 


0.06 ± 0.01 
0.10 ±0.03 


averaged 
values 


^ 


1 
3 


_ 


5.6 ± 0.4 
6.5 ±0.1 


0.7 ± 0.2 
0.36 ± 0.06 


2 
4 


~ 


7.2 ±0.1 
7.02 ± 0.07 


0.09 ± 0.02 
0.16 ±0.05 



For large and extended magnetic-field configurations, such 
as sunspots and active regions, this may not be a problem due to 
a slow er evolution of the magnetic field (see e.g. ISpadaro et al.l 
|2004 who state that the lifetime of active regions range from 
days to months), but for small scale magnetic fields (see e.g. 
a comparison between plage and quiet Sun, Title et al. 1992), 
which evolve in the range of minutes or even faster, this certainly 
is one. This problem can be dealt with in two ways. One way to 
solve this shortcoming is that one reduces the available FOV (as 
for data sets II, III & IV) to get a better correspondence between 
filtergram data and magnetogram data. A smaller FOV can be 
scanned faster, and also the cadence of the spectro-polaiimetric 
scans can be increased. On the other hand, a smaller FOV gives 
a statistically less robust result than a larger FOV would give. 
The other possibility would be to take several filtergram images 
of the whole FOV and take only a stripe of the image that is 
co-temporal with a certain stripe of the magnetogram. This pro- 
cedure would be suitable for data set I, but would cause a major 
(not practicable) eff'ort of data alignment since one should keep 
in mind that the stripes should be aligned on sub-pixel level, 
which gets more and more complicated, the smaller the FOV 
becomes (less alignment informatiorQ). 

Therefore we decided to take one co-aligned filtergram and 
compare it with the full magnetogram. This gives just a rough 
estimate of the true magnetic field strength distribution caused 
by the temporal shortcomings mentioned before and by the fact 
that the evolution of MBPs is fast and that they show a short life- 
time in the range of minutes (e.g. Ide Wiin et al.ll2005l lUtz et al.l 
[2010). On the other hand, the result for data set I was verified by 
the three other, independently taken data sets. These congruent 
results are gained despite the poor temporal coincidence between 
magnetic field information and filtergrams of data set I. This ob- 
vious antagonism leads us to the conjecture that the underlying 



' Less alignment information implies more noise and errors intro- 
duced in the analysis by increased alignment errors. 



magnetic field changes on a different (much larger) characteristic 
timescale than the MBP fea tures do. This can be also c oncluded 
by comparing two works of Ide Wiin et al.l (I2005ll2008h . who in- 
vestigated in one study the lifetime of MBP features deriving a 
mean lifetime of 3.5 min and in the other study the authors in- 
vestigated patches of magnetic fields and found longer lifetimes 
with a mean of about 10 min. These observational constraints 
give lise to the idea that MBPs are formed over pre-existing 
magnetic field patches. These field patches are more extended 
than the associated MBP that increases the probability of mea- 
suring a magnetic field strength signal even after a certain time 
lapse. After a typical lifetime in the range of minutes, MBPs 
dissolve and can reappear later on at the same spot. Whether 
this reoccuiTence really happens has to be investigated in more 
detail in future works. For data set I (sunspot group), we expect 
strong magnetic fields and large magnetic fluxes (sunspot fields), 
thus the underlying magnetic field may be even more stable (in 
the range of an hour), and several recurrent MBPs might be ex- 
pected. This may explain the good agreement of the results for 
the different data sets. 



Probability distributions: In this paper we employed a quite 
simple probability distribution with six free parameters consist- 
ing of two log-normal distributions in a first step. As a result the 
derived parameters are two geometrical mean values /u, two mul- 
tiplicative standard deviations cr, and two weighting coefficients 
w. As one can see by the corresponding distributions (Figs.[3][5] 
[7] and[8]l, this simple model does not fit the observed distribution 
perfectly, which can also be concluded from the corresponding 
reduced ;t'^ values, which show values between 3.0 and 8.3. The 
partly high x^ values can be explained e.g. by the impossibil- 
ity of such an easy fitting model to differentiate and fit the two 
distinct bumps at about 1000 G and 1500 G, respectively. In the 
following we discuss more complicated fitting models by adding 
first another (third) log-normal component (to enhance the fitting 
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for intermediate fields) and finally a fourth log-normal compo- 
nent (splitting the collapsed state into a weak and a strong field 
case). Furthermore we want to offer some possible explanation 
for the appearance of these components. 

Starting with data set I and applying the medium-complex 
fitting model on the results (three log-normal components) 
shown in the third panel from top of Fig.|9] we see that the agree- 
ment between the measured distribution and the fitted curve is al- 
ready better when compared to the simplistic model (given in the 
second panel from top). The solid black line shows the resulting 
total fit, while the single components are shown in color (red for 
the backrgound, blue for the pre-collapse or intermediate state, 
and purple for the collapsed MBP state). The quality of the fit 
improved from a value of 3.9 to 1.3. The figure also shows that 
the fitted distribution still misses a major feature around 1 000 G. 
To accomplish a well-fitted model distribution, we finally added 
a fourth component. The plot of this four-component model is 
shown in the bottom panel. One can see that the four-component 
model fits the measured distribution nearly perfectly. This is also 
expressed by the fit-quality measure, namely the reduced ;i'^ that 
now has a value of 0.8. 

Now the question arises whether we can assign an intuitive 
meaning to the three- and four-component models, respectively. 
We tried to do so and came up with the following interpretations 
for the components (see also the legend in the figure): 

- background: represents probably uncollapsed fields and 
purely background magnetic fields related to the acting so- 
lar dynamo processes; shows a geometric mean value of 
(270 ■ I H- 2) G. 

- pre-collapse: consists of MBPs in their pre- and post- 
collapse phase (an intermediate stage in which background 
magnetic fields get unstable and finally reach the collapsed 
state and a new equilibrium; obtained geometric mean value 
of (670 -K 1.4) G). 

- weak collapsed fields: the collapsed state of flux tubes. We 
assume that the weak component refers to collapsed intranet- 
work fields (estimated geometric mean value of (1120 ■ 
I- 1-2) G). 

- strong collapsed fields: in the case of the four-component 
model, the collapsed state seems to split up in a weak and 
a strong component. The strong component probably refers 
to collapsed network magnetic fields featuring a geometric 
mean value of (1340 -In- 1.1) G. 

From this interpretation it is clear that the contribution of the 
components in each single data set may vary depending on 
the kind of observed solar surface (e.g. exhibiting more or 
less magnetic network). Nevertheless, if the suggested four- 
component model is indeed related to true physical processes, 
it should be possible to fit all four data sets with the same model. 
Furthermore, the fitted distribution functions should yield nearly 
the same parameters, and only the weighting of the components 
should change from case to case. To test this we applied the 
medium complex and sophisticated model also to all of the other 
data sets (shown in Figs. [10] [TT] and[T2li. 

It is indeed possible to fit all four data sets very well with 
the three- and four-component models. The quality of the fit in- 
creases (see the figure captions), as can be supposed due to the 
higher number of fitting parameters. Only in the case of data set 
II also the simplistic fit gives a good agreement with the mea- 
surements expressed by reduced x^ of 3.0. Nevertheless also in 
this case the medium complex and the sophisticated model fit 
better by far. An overview of these parameters can be found in 
Tables |6] and |7] One can see that the parameters describing the 



different components are obtained in all cases with very similiar 
values. This is an indication that the distribution actually con- 
sists of several components and that the easier model misses the 
weaker componento 

We can conclude that the two component model akeady 
reasonably fits the observed data set, and yet a three- or four- 
component model fits the observed distributions significantly 
better The necessity of the additional components can be in- 
directly verified by the small error in the weighting parameters 
(see Tables |6] and |7]). In case of overfitting or introducing unnec- 
essary components to a fit, the error in the weighting parame- 
ter should be in the range of the weight itself. Furthermore, we 
can assign by intuition a distinct meaning to each component of 
the complex models. Nevertheless, whether such a model makes 
sense would have to be proven by theoretical means and by in- 
vestigations of single MBPs in more detail; i.e., do they form out 
of a background magnetic field and go through an intermediate 
state to finally reach the collapsed state? What happens there- 
after, and can they reach two different distinguishable collapsed 
states (associated to network and intranetwork fields)? 

In the future more detailed investigations will either prove 
the physical significance of these complex models or show that 
they are rather complicated functions that fit the observed dis- 
tributions without any physical meaning buried in their compo- 
nents. 

Finally, we stress that our sophisticated fit model is highly 
non-linear, and thus small changes in the starting conditions can 
change the outcome for the components (especially the weak 
ones). It is therefore necessary to give an educated guess about 
the right starting parameters. This means, e.g., visually identify- 
ing the peaks in the distribution and setting the starting parame- 
ters for the geometrical mean value of the distributions accord- 
ingly, i.e the value should be centred on the peak. Nevertheless, 
the main result (the bump around 1300 G in the magnetic field 
distribution for MBPs) can be deduced from the sophisticated 
models, as well as from the easier model, and has to be taken as 
evident. 



6. Conclusions 

We analysed four different data sets of Hinode/SOT compris- 
ing G-band and SP data. The four data sets show different so- 
lar surface features (active region with sunspots and pores, and 
quiet Sun). The magnetic field strength distribution of MBPs 
was obtained for those four data sets by comparing fully in- 
verted SP data (magnetograms) with carefully co-aligned G- 
band filtergrams. All the data sets showed the same magnetic 
field strength distribution behavior with a strong peak around 
1300 G. This value was theoretically predicted by the convec- 
tive collapse model (as e.g. stated by Spruit 1976) and is now 
observationally verified in this paper Furthermore, we showed 
that the magnetic field background signal (non-collapsed flux 
tubes), as well as the complete distribution, can be fitted by log- 
normal distribution components. The number of components for 
the model and the interpretation are outstanding, but probably 
linked to the dynamo and formation processes that form the mag- 
netic fields/MBPs. The log normal component related to the col- 
lapsed flux tubes shows a contribution ranging from 10% to 30% 
(except for the active region case where it reaches 60%), indicat- 



'" If components would not be required or randomly fitted into the 
distributions, they should not be assigned similiar values for each data 
set but would instead turn out to be assigned difi^erent values for each 
case. 
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ing that a large fraction of MBPs exhibit weak and uncollapsed 
magnetic fields. 
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